Monitoring stress related velocity variation in concrete with a 2.10 
relative resolution using diffuse ultrasound. 
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Ultrasonic waves propagating in solids have stress-dependent velocities. The relation between stress 
(or strain) and velocity forms the basis of non-linear acoustics. In homogeneous solids, conventional 
time-of-fiight techniques have measured this dependence with spectacular precision. In heteroge- 
neous media like concrete, the direct (ballistic) wave around 500 kHz is strongly attenuated and 
conventional techniques are less efficient. In this manuscript, the effect of weak stress changes on 
the late arrivals constituting the acoustic diffuse coda is tracked. A resolution of 2.10 -5 in relative 
velocity change is attained which corresponds to a sensitivity to stress change of better than 50 
kPa. Therefore the technique described here provides an original way to measure the non-linear 
parameter with stress variations on the order of tens of kPa. 



PACS numbers: 43.20.Jr, 43.25.Ed, 43.25.Dc, 43.20.Gp 



Concrete is a complex heterogeneous mix of various 
ingredients of different sizes ranging from centimetric to 
millimetric gravel and sand, to the micrometric porosity 
of the cement paste. Over the last thirty years or so, 
various sonic and ultrasonic techniques have been devel- 
oped to image, characterize or control the materiali^^. 
Applications mostly work in the low frequency regime (a 
few Hz up to 50 kHz), where the wavelength is greater 
than the sizes of heterogeneities. In this regime, direct 
waves do not suffer too much from scattering by the het- 
erogeneities and wave arrivals can be simply related to 
ballistic (direct or simply reflected) propagation paths. 
Nevertheless, because the wavelengths are much greater 
than the microstructure of the concrete, the fine details 
like micro-damage or millimetric cracks are hardly visi- 
ble. In the high frequency regime (above 50 kHz), the 
elastic wavefield is sensitive to these small details, but 
the ultrasonic wavelength is also of the order of the ag- 
gregates' size. This results in a strong attenuation of 
the direct (ballistic) wave££, and the appearance of a 
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FIG. 1. Schematic view of the experimental setup. 



long lasting coda made of waves that have been multiply 
scattered by the internal microstructure^. The multi- 
ple scattering regime is often associated with a catas- 
trophic loss of information, such that most conventional 
imaging or monitoring techniques fail. Nevertheless, as 
diffuse waves travel much longer paths than direct or 
simply reflected ones, they are much more sensitive to 
weak perturbation of the medium. This idea has been 
exploited in geophysics to quantify weak changes in the 
velocity of the earth crusl£: by comparing two seismic 
coda obtained from fixed source and receiver at two dif- 
ferent times, it is possible to monitor weak velocity vari- 
ations in the medium. This technique is known as the 
seismic doublet technique or (more recently) coda wave 
interferometrylQ. 

Various phenomena can cause acoustic velocity vari- 
ations in concrete, including chemical reactions, and 
changes in water content, temperature or stress. In this 
work we focus on the latter effect, i.e., the relation be- 
tween stress/strain and acoustic velocity, which is known 
as the acousto-elastic effec t 11 ' 12 ! 13 . Under vertical uni- 
axial stress Act the ultrasonic velocities in a sample will 
evolve as (to first order): 



da 



Aa + o(e), 



(1) 



where ij stands for the wave velocity in direction i for 
a particle motion in j. The dimensionless non-linear pa- 
rameter j3 is defined as: 
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and depends on the second (Lame) and third (Mur- 
naghan) order elastic constants. E w 3.9.10 10 Pa is the 
Young's modulus. There has been much less attention 
on the c-dependence of V in heterogeneous media, 
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FIG. 2. Typical waveform h[t] collected through the concrete FIG. 3. Relative velocity variations dV/V versus stress vari- 
sample. The processed time-window marked by two vertical ation Act. 
dashed lines extends from 50 /is to 500 /is. 



such as concrete---, than in homogeneous media, and 
most results have considered variations under significant 
applied stress. To our knowledge, the use of diffuse 
ultrasound to evaluate the effect of stress loading on 
velocity variations has only been reported by Gret 
et alJ£ and Payani^. In these works the multiple 
time-window doublet technique was used to infer relative 
velocity variation with a relative resolution of 10~ 3 to 
10~ 4 for stress variation of the order of 1 MPa. Here, we 
propose to monitor weaker relative velocity variations in 
concrete down to 2.10 -5 under weaker load variations 
(50 kPa). Such precision requires high frequency data 
(here 500 kHz diffuse ultrasonic data) and a more 
sensitive data processing technique than the doublet 
technique. 

The experimental data presented here were acquired on 
a concrete sample of 16 cm in diameter, and 28.5 cm in 
height. The concrete was prepared with 17% (in weight) 
of cement, 31% of fine sand, 43% of gravel (size rang- 
ing from 2 to 10 mm) and 9% of water. Two ultrasonic 
transducers are used in the experiment, one as a source 
and one as a receiver (labeled S and R, respectively, in 
Fig. [T|) . The lateral size of the active part of the trans- 
ducers (0.7 mm) is much smaller than the wavelength, 
which makes them very sensitive to multiply scattered 
waves. The transducers are glued onto the sample using 
a hot chemical glue (phenyl-salicylic acid) that solidifies 
with cooling (below 43°C). The ultrasonic experiment 
was performed several months after casting the sample. 

During the ultrasonic measurements the concrete spec- 
imen was subjected to uniaxial loading. A preliminary 
load was applied (5 MPa) at which point the reference 
ultrasonic data ro(t) were acquired. The applied uniax- 
ial load was then increased to 5.5 MPa in ten step of 
50 kPa with ultrasonic acquisition at each stress incre- 
ments, k, to provide rk(t). Each ultrasonic acquisition 
involved a source excited with a 10V chirp s(t) of fre- 
quencies ranging linearly from 200 kHz to 800 kHz and 
of duration 10 s. The received signal rk(t) was ampli- 
fied by 40 dB and stored in the computer. The recorded 
waveforms were subsequently cross-correlated with s(t) 
(time-compression) to provide an estimate of the impulse 



response of the sample hk(t) = Tkit) x s(t), for each incre- 
ment k. The effective duration of the time-compressed 
signal hk(t) is about 0.6 ms. ho(t) is the reference wave- 
form corresponding to r${t) correlated with s(t). It is 
noted that the source and receiver positions were kept 
fixed for all acquisitions and the excitation waveform s(t) 
was not changed. Note also that the velocity change due 
to strain is one order of magnitude smaller than the one 
due to stress change, and is neglected here. The time 
separating two consecutive measurements is 5 minutes, 
over which no stress change was recorded and the ultra- 
sonic records were reproducible. Creep was observed over 
much longer time scales (several hours) and was negligi- 
ble in the present experiment. 

A typical ultrasonic record is plotted in Fig. [51 This 
record is composed of: (1) the impulse response of the 
concrete sample; (2) a peak at t=0 marking the initial 
excitation and due to the electronic cross-talk between S 
and R; (3) background electronic noise (about 1% rela- 
tive amplitude). The impulse response consists of a di- 
rect wave arriving at 0.023 ms, which is barely visible, 
and a following coda due to strong multiply scattered 
waves. The average intensity of the coda, I(r,t), can be 
described by a diffusion envelope: 



I(r,t) 
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(3) 



where Iq is the energy released by the source, D is the dif- 
fusion constant, r is the source-receiver distance and 7 is 
the dissipation rate. In this experiment around 500 kHz, 
we find that y/l best fits the data for D = 5 mm 2 / us 
and 7 = 25 ms -1 . This will represent a slight under- 
estimate for D because Eq. [3] only holds for unbounded 
media. Nevertheless, D and a values are in agreement 
with the literature^ 8 ., and our rough estimation of D 
confirms that we are in the strong multiple scattering 
regime. 

For each load k, the record hk(t) is compared to the 
reference waveform ho (t) to evaluate the relative velocity 
change in the concrete sample. Because of strong mul- 
tiple scattering, the P- and S-wave arrivals mix together 
and rapidly attain the equipartition regime^. Let ay- 
be the relative fraction of time that the wave spends in 
the state ij (polarization in the direction i while prop- 
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agating in the direction j), then the velocity variation 



is dV — (aij^^-jlS.a. Note that we have J2 a v = 1 

and that (.) denotes ensemble averaging. If the waves 
are statistically equipartitioned in all directions, we have 
a,i=j — 9% and ^ a^j = 91%—. In general, ay has a 
more complex distribution, but in this work we focus on 
the average quantity (a^/Jy). For simplicity, we also as- 
sume that the waveforms are just stretched in time. This 
is again a first order approximation that applies to ex- 
perimental data in the considered time-window. Never- 
theless, velocity variations will progressively decorrelate 
the waveform^. This weak decorrelation, also called dis- 
tortion, is not studied here although it contains precious 
information on the medium and its evolution. In partic- 
ular, we stress that strong distortion and weak stretching 
might indicate a change of structure rather than a global 
velocity change. However here we do not see significant 
distortion. To summarize, the velocity change dV for 
load k results in: 



h k [t]=ho [t{l + v k )]+n(t) 



(4) 



with v k = dV/V and n(t) accounting for additional fluc- 
tuations, including electronic noise and the small distor- 
tion mentioned earlier. Two processing techniques have 
been proposed in the literature to estimate dV/V. The 
first one, called the seismic doublet technique, was devel- 
oped for geophysical purposes about twenty years agcA 
The idea is to measure a time-shift between two differ- 
ent records in a limited time-window centered at time 
t. By measuring such time-shifts (or delay) St at dif- 
ferent times t in the coda, it is possible to evaluate the 
velocity variation, which is simply the average slope of 
St(t): dV/V = —St/t. Such an approach implicitly as- 
sumes that the time-shift is constant within the consid- 
ered time-window, which is generally not the case. An- 
other ide a 18 i 19 is to interpolate h k at times f(l — v) with 
various relative velocity changes z^£. Therefore v k is the 
v that maximizes the cross-correlation coefficient: 
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Contrary to the doublet technique, we do not assume a 
constant time-shift in the considered time window [0 T]. 
The comparison between the two waveforms can there- 
fore be performed over the entire record at once, which 
results in a more stable, and thus more precise, esti- 
mation of dV/V. The comparison between the doublet 
and the stretching technique will be subject to further 
investigational. Assuming that ho and hi are stationary 
waveforms and are well described by Eq. [H we have a 
theoretical estimation of CC: 



CC k (y)=A- 



I A f P(f) sinc ( 27r / - Vk) T) df 
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with / the frequency, A/ the bandwidth, p(f) the spec- 
tral density. The constant A depends on the variance of 



h, noted (h 2 ) and the variance of the additional fluctua- 
tions, noted (n 2 ), 



A : 
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B(v) is a random process of zero mean and standard 
deviation: 



1 



(n 2 ) + 2^{n 2 ){h 2 ) 



(8) 



If the amplitud e of t he signal A is much greater than 
the fluctuations yj (B 2 ) of the cross-correlation coefficient 
CCfc, the maximum of CC k is obtained for the stretch- 
ing v — v k . It is interesting to note that the peak of the 
sine function is visible even if the fluctuations (or noise) 
n(t) are strong. In such a case, increasing the integra- 
tion time T or the frequency bandwidth A/ can compen- 
sate for this strong noise. This is a crucial advantage of 
the present technique compared to the doublet technique. 
The dV/V resulting from our measurements are plotted 
in Fig. [3] versus Act. A linear regression of the data yields 
(vifc) = 10~ 6 kPaT 1 , which gives an estimate of the 
non-linear parameter of our sample: [3 = (dij/3y) ~ —40. 
The standard deviation from the linear regression for the 
relative velocity change dV/V is 2.10~ 5 , which indicates 
that the precision on the evaluation of (3 is 5%. Note that 
during the whole experiment, the surrounding air tem- 
perature fluctuated within less than 0.8°C. Temperature 
variations are too weak and thermal diffusion too slow 
in the sample (compared to the duration of the experi- 
ment) to induce noticeable cil / 10 ) 18 i 22 i 23 . This assertion 
is confirmed by the fact that the velocity variations dV 
are totally decorrelated from the air temperature varia- 
tions. We therefore propose that the variations are due 
primarily to stress effects. 

A potential application of this technique is on-site as- 
sessment of the non-linear constant j3 of a concrete struc- 
ture. Our technique is sensitive enough to provide reli- 
able measurements with additional stress Act of the or- 
der of 100 kPa. As the non-linear parameter strongly 
depends on the micro-damage of the concrete, our tech- 
nique would form an efficient and non-destructive way to 
probe the integrity of civil engineering structures. An- 
other application is to evaluate the state of stress of a 
given structure with a predetermined (3 parameter, and 
can be extended to geophysical medi a 24 i 25 ' 26 . 
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